Ultralong-range polyatomic Rydberg molecules formed by a polar perturber 
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The internal electric field of a Rydberg atom electron can bind a polar molecule to form a giant 
ultralong-range stable polyatomic molecule. Such molecules not only share their properties with 
Rydberg atoms, they possess huge permanent electric dipole moments and in addition allow for 
coherent control of the polar molecule orientation. In this work, we include additional Rydberg 
manifolds which couple to the nearly degenerate set of Rydberg states employed in [S. T. Rittenhouse 
and H. R. Sadeghpour, Phys. Rev. Lett. 104, 243002 (2010)]. The coupling of a set of (n + 3)s 
Rydberg states with the n(l > 2) nearly degenerate Rydberg manifolds in alkali metal atoms leads 
to pronounced avoided crossings in the Born-Oppenheimer potentials. Ultimately, these avoided 
crossings enable the formation of the giant polyatomic Rydberg molecules with standard two-photon 
laser photoassociation techniques. 

PACS numbers: 33.80.Rv, 31.50.Df 



I. INTRODUCTION 

The physics of Rydberg atoms and molecules has devel- 
oped into a quasi-sustainable "ecosystem" over the last 
decade, mainly due to a) the exaggerated properties of 
Rydberg atoms and molecules (long lifetimes, large sizes, 
and scalability), and b) exquisite experimental control 
over these properties through the advent of ultracold 
atomic Rydberg samples [TJ [5] . Some landmark devel- 
opments over the last few years in the field of Rydberg 
physics have been the creation of a frozen Rydberg gas 
[3] , the demonstration of a Rydberg blockade scheme and 
the subsequent realization of a Rydberg qubit gate |U- 
7 , the formation of ultralong-range isotropic Rydberg 
molecule [8], the creation ultracold neutral plasmas [9], 
as well as the formation of highly-magnetized Rydberg 
antihydrogen atoms |10) . These rapid advances now al- 
low for few-body and many-body effects to be realized 
in laboratory settings and be prototyped for simulation 
of strongly interacting spin chains [TTHT5] . The other in- 
gredient of our proposed polyatomic Rydberg complex, 
namely, ultracold polar molecules, have also been touted 
as toy models for the simulation of many-body condensed 
systems and for the realization of quantum gates [T3]. A 
main appeal of polar molecules is that they may possess 
sizeable permanent dipole moments. 

In this work, we harness the long-range interaction of 
a Rydberg atom with a polar molecule (Rb is used as an 
ubiquitous example, but other alkali metal atoms would 
serve the purpose), to demonstrate that ultralong-range 
polyatomic molecules with enormous permanent dipole 
moments can form from combining a Rydberg atom and 
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a polar molecule. This system was first studied in [TS] 
where it was shown that the Rydberg electron can be 
used to "drag" the polar molecule into a preferred orien- 
tation and that this process can be controlled by a Ra- 
man microwave pulse scheme. Ultralong range Rydberg 
molecules in general came in vogue after a proposal that 
such molecules could form from zero-range interaction of 
electrons with ground state atoms [TB] . The recent exper- 
imental realization of a class of such molecules [S] has led 
to revived theoretical and experimental activity [T7J HH] . 



The interaction of the polar perturber with the Ryd- 
berg atom strongly couples the field-free atomic Rydberg 
states by means of the dipole's electric field. In [TS], the 
mixing of states was restricted to the nearly degenerate 
(negligible quantum defects) manifold of Rydberg states 
in Rb(n/ > 3). In the present work, we extend our previ- 
ous study by going beyond this degenerate perturbation 
theory approach. In this manner, the effect of additional 
Rydberg orbitals on the molecular Born-Oppenheimer 
(BO) potentials is probed. In particular, the proxim- 
ity of the n(l > 2) Rydberg manifold to a single (n + 3)s 
state introduces a strong coupling, allowing for the exper- 
imental preparation of the proposed polyatomic Rydberg 
molecules via a standard two-photon laser excitation. 



The paper is outlined as follows. In section |TT] the 
model Hamiltonian describing the polar molecule per- 
turber is introduced. The adiabatic Hamiltonian for the 
full Rydberg atom plus polar molecule complex is subse- 
quently provided in section |III| The resulting adiabatic 
potential surfaces are discussed in section |TV| In section 
V we discuss the admixture of s-wave electron character 
in the giant Rydberg molecule. We conclude by a brief 
summary and an outlook on further research directions 
in section IVI1 
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II. POLAR MOLECULE MODEL 
HAMILTONIAN 

Before we discuss the properties of the proposed poly- 
atomic Rydberg molecules, we must introduce an effec- 
tive Hamiltonian for the molecular perturber. In this 
context, we model the polar molecule as a two-level 
molecule in which the opposite parity states are mixed 
in the presence of an external field, i.e., 
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where A is the zero field splitting between the two molec- 
ular states (as in a A-doublet molecule), d is the perma- 
nent dipole moment of the molecule in the body fixed 
frame, and F ext is an electric field external to the polar 
molecule. Specifically, the external field stems from the 
Rydberg electron and the Rydberg core, i.e., 
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where e is the electron charge, R is the core-polar 
molecule separation vector and r is the position of the Ry- 
dberg electron with respect to the core. In other words, 
describes the coupling between the internal state of 
the polar molecule and the Rydberg atom. 

It should be noted that only depends on the magni- 
tude of the external field, which is appropriate for static, 
homogeneous fields. For a rigid rotor molecule, such 
as KRb, the rotation of the molecule happens on much 
slower time scales than the Rydberg electron orbital time. 
In this case, the coupling between the Rydberg atom 
and the polar molecule accounts for the rotation of the 
molecule, 
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where d is the rigid rotor dipole moment. This interac- 
tion leads to BO potential surfaces which depend on both 
R and the orientation of d with respect to R; this more 
complete treatment is the subject of an ongoing study. 
In contrast, a A doublet molecule, such as OH, has a per- 
manent dipole moment that arises from the interaction 
of two opposite parity electronic (e, /) states. Rotational 
transition energies in A-doublct molecules are usually or- 
ders of magnitude larger than the typical doublet split- 
ting energies. For the fields provided by highly excited 
Rydberg atoms, of order F ext ~ 10 -6 a.u., our model 
Hamiltonian thus provides an excellent description of 
such molecules. 

The electron-dipole interaction in and in the off- 
diagonal elements of has a critical value. If the dipole 
moment, do, is larger than the Fermi- Teller critical value, 
d c = 1.63 D, an infinite number of bound states form |19j . 
Furthermore, an electron scattering off a super critical 



dipole (d > d c ) is sensitive to the detailed, short-range 
structure of the polar molecule. To avoid the complica- 
tion of super-critical dipole scattering and electron trans- 
fer from the Rydberg atom to the polar molecule, we will 
deal exclusively with polar molecules whose dipole mo- 
ments are subcritical, i.e., do < 1.63 D. 



III. THE ADIABATIC HAMILTONIAN 

For simplicity, we assume that the external electric 
field is along the intermolecular axis, such that be- 
comes 
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where 



- Rj ■ R/R r~R. This means that 
the projection m of the Rydberg electron angular mo- 
mentum along R is conserved. To find the BO potentials, 
we solve the adiabatic Schrodinger equation at fixed po- 
lar molecule location R = Rz, 



H ad ^ {R; r, a) = U(R)^(R; r, a), 



with 



H a d = Ha + H, 
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The first term, Ha 



i a = ~2^V^ + Vi(r), describes the 
unperturbed Rydberg atom; the core penetration, scat- 
tering, and polarization effects of its valence electron are 
accounted for by the /-dependent model potential Vi(r) 
|20) , giving rise to the quantum defects of the low angu- 
lar momentum Rydberg states. m e is the electron mass, 
-0 is the electron wave function, and a is a coordinate 
signifying the internal states of the polar molecule. The 
eigenvalues U (R) serve as the sought-after BO potentials 
for the polar perturber. To solve 0, the total wave 
function is expanded in the basis {tpnim (r) |±)} where 
(r) is an unperturbed Rydberg orbital and |±) are 
the polar molecule parity states. The atomic and molec- 
ular degrees of freedom are coupled by the electric field 
that mixes both the Rydberg orbitals as well as the 
parity states of the polar molecule, cf. 0. 

In [T5], the BO potentials were found using degener- 
ate perturbation theory in the (nearly) degenerate set of 
Rydberg orbitals {ip n (l>2)o (f)}- The molecular Hamilto- 
nian Q can be prediagonalized to yield the eigenvalues 
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where F c — A/2rfo is the critical external field strength 
at which the molecule becomes completely polarized. In- 
serting in for F ext the i?-dependent eigenstates of the 
Rydberg atom exposed to the field results in the po- 
tentials found in |15j . A more complete treatment of the 
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system requires that the total adiabatic Hamiltonian is 
diagonalized in a complete set of Rydberg orbitals. Be- 
cause in the latter case, the orbitals are in general not 
degenerate, the prediagonalization scheme cannot be di- 
rectly used beyond the degenerate perturbation theory. 

In the present work, we go beyond the degenerate per- 
turbation theory framework adapted in |15| . To this end, 
the Hamiltonian ^ is diagonalized in an extended basis 
set comprising several Rydberg n manifolds. By con- 
struction, the atomic Hamiltonian Ha is diagonal in the 
Rydberg basis {^ (r)}, 



The off-diagonal matrix elements of ^ are due to H mo i 
which couples the Rydberg electron states to the internal 
parity states of the polar molecule. They are determined 
by evaluating the integrals 
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An intuitive closed form of (|Sj) is obtained if we expand 
the electron electric field contribution in spherical har- 
monics. Using the multipole expansion of [2"Tj yields 
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where fii is the quantum defect for the l-th partial wave 
Rydberg state and a denotes the parity state of the polar 
molecule. The quantum defects used in this work are 
those for the rubidium atom [2]: fi s — 3.13, fi p = 2.65, 
Hd = 1-35, (if = 0.016, and ^/ >3 s» 0. While our focus is 
on Rb, the method can be extended to any highly excited 
Rydberg atom by using the appropriate quantum defects. 
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We note that, with some generalization, this procedure 
can also be used to evaluate the interaction matrix ele- 
ments between a permanent dipole and the Rydberg elec- 
tron, cf. (J3j) . Using ([£]), the off diagonal matrix elements 
of (§ are given by 
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where R n i (r) is the radial wave function for an electron 
in the nl Rydberg state. 

The adiabatic Hamiltonian matrix assumes the follow- 
ing form, 
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where I is the unity matrix, ERyd is the diagonal matrix 
representation of Ha, and F is the electric field matrix 



whose elements are given by ( 10 ) with a diagonal offset of 
1 /R 2 due to the contribution of the electric field from the 
positive ionic core to the total electric field, as defined in 
Q. The BO potentials are found by diagonalizing the 
matrix (11) at each R. Special attention must thereby 



be drawn to the actual size of the basis set in order to 
ensure convergence of the potentials. 



IV. BORN-OPPENHEIMER POTENTIALS 

In [T5], it was assumed that for I > 2, the Rydberg or- 
bitals of Rb are degenerate and that only this degenerate 
set of states is required to converge the BO potentials. 
Here we extend this treatment to including additional 
Rydberg orbitals and the small, but finite, /-wave quan- 
tum defect. Figure [l] compares the BO potentials using 
degenerate perturbation theory (solid curves, as in |15j ) 
to those found by including the /-wave quantum defect 
(dotted curves) for the n = 25 state of rubidium; a polar 
molecule dipole moment of do = 0.40 a.u., and a zero- 
field splitting of A = 1.85 x 10~ 7 a.u. is considered. 

We reiterate some of the most salient features of the 
proposed molecules. The modulations that form a se- 
ries of wells reflect the oscillatory nature of the Ryd- 
berg electron wave function. The outer most wells in 
the lowest two potentials are deep enough to support 
many vibrational levels. The resulting giant polyatomic 
Rydberg molecules share several features with previ- 
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FIG. 1: The BO potentials for the Rb(n = 25) Rydberg atom 
and a polar molecule are shown. The potentials are calcu- 
lated at the level of degenerate perturbation theory including 
(dashed red curves) and ignoring (solid black curves) the /- 
wave quantum defect of the Rydberg electron. The polar 
perturber in this example is a molecule with a dipole moment 
do = 0.40 a.u. and a zero field splitting A = 1.85 x 10~ 7 a.u. 



ously predicted homonuclear Rydberg molecules, the so 
called "trilobite" molecules [111 [221 123]- The size of 
these molecules scales as R ry d oc n 2 and the well depths 
scale as Vp cx do/n 3 . Unlike for homonuclear Rydberg 
molecules, the anisotropic nature of the electron-dipole 
interaction creates two different internal configurations. 
Corresponding to the well labelled R in figure [I] in one 
configuration the dipole of the polar molecule is oriented 
towards the positive core. In the other configuration - 
labelled L - the dipole of the polar molecule is oriented 
away from the positive core. 

As anticipated, inclusion of the small /-wave quantum 
defect only modifies the potentials close to the dissocia- 
tion threshold limit, i.e., near the Rb(n = 25) limit. In 
the lowest wells, where the Rydberg molecules form, the 
two sets of potentials are slightly shifted with respect to 
each other while being otherwise almost identical, indi- 
cating that including the /-wave quantum defect has al- 
most no influence on the behaviour of the resulting giant 
molecule, nor on the polyatomic molecular dipole mo- 
ment. 

By including more Rydberg orbitals in our basis set, 
we can explore the range of validity of the degenerate 
perturbation theory approach. For smaller dipole mo- 
ments, do < 0.4 a.u., relatively few basis states are re- 
quired to achieve convergence. In fact, the converged 
potentials are only slightly different from those shown 
in figure [T] As the dipole moment increases, more and 
more Rydberg orbitals arc required. This is due to the 
localizing effects of the dipole-electron potential. For a 
nearly critical dipole moment, do ~ 0.63 a.u., the elec- 
tron is almost entirely localized at the location of the 
dipole. To accurately describe this (small angle) local- 



FIG. 2: (a) The BO potentials for the Rydberg-polar 
molecule system are shown, calculated including lower angu- 
lar momentum states using a polar molecule A-doublet parity 
splitting of A = 1.85 x 10~ 7 a.u. and a dipole moment of 
do — 0.40 a.u. (subcritical dipole). (b) The same as (a) with 
a dipole moment of do = 0.57 a.u. (the permanent dipole mo- 
ment of CD). The zero in energy is the Rb(n = 25) threshold. 
The 28s level and zero field molecular splitting have also been 
labelled. 

ization requires an inordinately large number of Rydberg 
orbitals. Furthermore, for larger values of do, the elec- 
tron wave function penetrates into the short-range re- 
gion of the perturber molecule and probes the detailed 
structure of the molecule. We therefore restrict here to 
sub-critical dipole moments, do < 0.6 a.u. Because the 
Rydberg spacing scales as n~ 3 , larger sets of orbitals are 
required to converge the potentials attached to higher 
Rydberg thresholds. The calculations reported here are 
for Rb(n = 25) manifold of states, though the qualitative 
behaviours we discuss will persist for higher n. To acquire 
converged potentials for n = 25 at the largest dipole mo- 
ments considered here, do = 0.60 a.u., our basis includes 
all electron angular momenta for n = 19, 20, 30, 31 as 
well as the 32(s,p, d), 33(s,p) and 34s Rydberg orbits, 
yielding a total of 331 electron basis states. 

Figure [2] shows the converged potentials attached to 
the n — 25 Rydberg level. As in figure [T] the polar 
molecule in figure [2]ja) has a dipole moment do = 0.40 
a.u. and a zero field splitting A = 1.85 x 10~ 7 a.u. (A 
has been chosen to be the A-doublet splitting of CD). 
Comparing the two figures shows that the molecular po- 
tentials are unaffected by the presence of the non-binding 
s-wave states for cases when the molecular dipole is less 
than ~ 0.4 a.u. The two curves converging to the Rb(28s) 
threshold correspond to the two opposite parity states of 
the polar molecule and are correspondingly split by the 
zero field splitting A. 

For larger dipole moments, it is essential to include 
the lower angular momentum states in our treatment. 
While for an unperturbed Rydberg atom, the latter are 
isolated from the degenerate I > 2 manifold, increasing 
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the dipole moment has the interesting effect of forcing 
the potential wells of the n = 25 to cross through the 28s 
potentials. This is possible because the Rb(ns) quantum 
defect (/x s = 3.13) has a small non-integer fraction, plac- 
ing the Rb(28s) close to the n = 25 degenerate manifold 
of states. In figure [2jb) , a corresponding example with 
d = 0.57 a.u.= 1.46 D - the permanent dipole moment 
of CD - is provided. In this case, the interaction be- 
tween the n = 25 molecular states and the Rb(28s) state 
is significant. This is generally true of all Rb(n(7 > 2)) 
degenerate manifolds and (n + 3)s interacting states. As 
a result and due to strong avoided crossings, potential 
wells capable of supporting bound vibrational levels ap- 
pear in the BO curves attached to the Rb(28s) threshold. 
Figure [3^ a) provides a close-up of the four BO potentials 
directly involved in the binding of the polyatomic Ryd- 
berg molecule. The vibrational wave functions for two 
such bound states in the outermost wells are also shown. 



V. S-WAVE ADMIXTURE 

The BO potential curves in figure 1 (also in [T5]) in- 
clude a nearly degenerate superposition of atomic or- 
bitals. As such, they solely contain contributions from 
high angular momentum states (I > 2) and hence are 
only accessible in experiments which " photoassociate" 
these molecules, if and when low-lying angular momenta 
are admixed into the degenerate manifold. The strong 
interaction of the (n + 3)s Rydberg state with the lowest 
BO potential curves belonging to the n(l > 2) manifolds 
in figure [3^b) admixes large amount of s-wave character 
into the wave functions. 

At the potential minimum, the adiabatic channel func- 
tion is well approximated by 
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i> ( r , X) = a (d ) ip d (r) \%d) + b (d ) tp s (r) \x s 
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where tpd and \xd) are the electronic wave function and 
internal polar molecule state, respectively, that include 
higher electron angular momentum with I > 2. ip s and 
\Xs) are the s-wave Rydberg electron wave function and 
the corresponding molecular state, respectively. The ex- 
pansion coefficients a (do) and b (do) depend on the dipole 
moment of the polar perturber as well as on its position 
R. Figure [3jb) gives the s-wave electron contribution, 
|6(<io)| 2 , as a function of the intermolecular separation 
distance for each of the adiabatic channel functions cor- 
responding to the potentials shown in figure [3ja). From 
figure [3^b) it can be seen that near the minima of the 
outer wells, the electronic state of the giant molecule has 
a significant s-wave character, approximately 30 — 40%. 
This implies the possibility that these molecules could be 
formed in a simple two photon excitation scheme similar 
to that used in the experimental realization of homonu- 
clear Rydberg molecules [8]. The n = 25 wells shown 
in figure [3]ja) are deep enough to support bound vibra- 
tional states, the radial wave functions of the first two 
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FIG. 3: (a) The four BO potentials, directly involved in 
molecular binding, are shown for the same parameters as in 
figure[2|b). Two vibrational wave functions are also shown at 
their binding energies in each of the outermost wells, (b) The 
s-wave character of the molecular wave function is shown for 
each of the four potentials. 



such states are shown at their corresponding binding en- 
ergies in the outer two wells. These molecular levels are 
well isolated and have readily accessible vibrational fre- 
quencies. 

The added s-wave character in the electron wave func- 
tion decreases the charge localization exhibited in the Ry- 
dberg molecule slightly. Even with this decrease, these 
Rydberg molecules exhibit massive dipole moments. Us- 
ing the scaling behaviour of [TS], the Rydberg molecule 
dipole moment d^yd can be found, 

dRyd ~ 1-3 [a (d )} 2 n 2 . 

For the case shown in figure [3]ja) this yields dn y d ~ 
1400D, a truly large dipole moment. The presence of 
these enormous dipole moments indicates that the poly- 
atomic Rydberg molecules could be sensitively controlled 
through the use of small external electric fields. Due to 
the extreme sensitivity of Rydberg electron to external 
fields, the behaviour of the BO potentials under the in- 
fluence of such fields is not immediately obvious and is 
the subject of ongoing studies. 

In figure |4^a-c), we provide snapshots of additional 
BO potentials responsible for molecular binding, namely, 
for d = 0.30, 0.55, and 0.60 a.u., respectively. As the 
dipole moment of the polar perturber increases, the L 
and R wells are pulled down through the s-wave Ryd- 
berg threshold. At do = 0.30 a.u., the well structure 
is not perturbed at all by the s-wave potential. For 
d = 0.55, the potential wells are strongly distorted by 
the presence of the lower threshold, while for do = 0.60 
a.u. the outermost well structures have passed through 
the s-wave threshold and correspondingly have a much 
smaller s-wave contribution. Figures Qd-f) show the s- 
wave electron contribution for each potential. For 0.5 
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FIG. 4: The variations of the BO potential curves with 
the strength of the polar molecule dipole moment do- ( a ) 
do = 0.30 a.u., (b) do = 0.55 a.u., and (c) do = 0.60 a.u. The 
respective s-wave contributions to the molecular wave func- 
tions are shown in panels (d-f). A A-doublet parity splitting 
of A = 1.85 x 10 -7 a.u. is considered. 



a.u.< g?o < 0.6 a.u. the s-wave character at the position 
of the outermost wells varies approximately 10 — 40%, 
enough to give a fairly large Rabi frequency for the cre- 
ation of Rydberg molecules in a standard two photon 
excitation scheme. 



interaction creates a series of BO potentials with an oscil- 
lating series of wells that reflect the Rydberg oscillations 
in the electron wave function. By extending the work 
in 15] to beyond the degenerate perturbation theory, we 
have shown that the small, but finite /-wave quantum 
defect changes little the behaviour of the resulting giant 
molecules. 

By including lower electron-angular momentum states, 
a new set of potentials attached to the (n + 3) s Rydberg 
level appear. These added potentials have no significant 
influence on electronic state of the Rydberg molecule for 
smaller subcritical dipole moments, do < 0.4 a.u. For 
nearly critical dipole moments, do 0.6 a.u., a series of 
avoided and level crossings is formed between the poten- 
tials formed from the degenerate n(l > 2) manifold and 
the (n+3) s-wave potentials. These avoided crossings lend 
significant s-wave character to the Rydberg electron state 
opening the possibility of creating the giant polyatomic 
Rydberg molecules using standard two-photon Rydberg 
excitation schemes. The series of avoided and level cross- 
ings create a complex structure of couplings between the 
various BO potentials. Extensions of this work beyond 
a simple two-state polar molecule to incorporate rigid 
rotor- type polar molecules, such as KRb, poses a sig- 
nificant challenge. The added angular behaviour of the 
polar molecule will create a set of two-dimensional po- 
tential surfaces which couple the rotational behaviour of 
the polar molecule to the vibrational state of the giant 
Rydberg molecule. Creating and examining this intricate 
energy landscape will be the subject of future work. 



VI. CONCLUSIONS AND OUTLOOK 

In this paper we have examined the formation of giant, 
polyatomic Rydberg molecules consisting of a A-doublct 
polar molecule and a Rydberg atom. We provided a more 
complete description of the polar molecule model initially 
used to predict these nearly macroscopic molecules [To] . 
By using a two-state model that is sensitive to external 
fields, a Rydberg electron is coupled to the internal state 
of the polar molecule. The resulting electron-molecule 
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